Electronic structure and electric-field gradients analysis in Celn^ 
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Electric field gradients (EFG's) were calculated for the Celn^ compound at both ^^^In and 
^^^Ce sites. The calculations were performed within the density functional theory (DFT) using 
the augmented plane waves plus local orbital (APW+lo) method employing the so-called LDA+U 
scheme. The Celns compound were treated as nonmagnetic, ferromagnetic, and antiferromagnetic 
cases. Our result shows that the calculated EFG's are dominated at the ^''"Ce site by the Ce-4f 
states. An approximately linear relation is intuited between the main component of the EFG's 
and total density of states (DOS) at Fermi level. The EFG's from our LDA+U calculations are 
in better agreement with experiment than previous EFG results, where appropriate correlations 
had not been taken into account among 4f-electrons. Our result indicates that correlations among 
4f-electrons play an important role in this compound and must be taken into account. 

PACS numbers: 71.20.-b, 71.28.+d, 76.80.+y, 71.27.+a, 71.20.Eh, 75.30.Mb, 75.25.+Z, 31.30.Gs, Sl.S.Ar 



I. INTRODUCTION 

Hyperfine interactions provide sensitive physical quan- 
tities such as electric field gradients (EFG's), which 
can be used to shed light experimentalljii'^ and 
theoretically^i'^ into the electronic states of materials. In 
this paper we have focused on the Celn^ as an interested 
system composed of strongly correlated 4f-electrons. It 
is a cubic heavy-fermion (HF) local moment antiferro- 
magnetic (LMAF) system at ambient pressure with a 
Neel temperature of 10.1 K— . This concentrated Kondo 
compound exhibits^ various fascinating and unexpected 
physical properties. The various properties of this heavy 
fermion originate from the fact that one cannot assign a 
definite localization to the 4f-states, irrespective of the 
applied conditions to the compound. The unexpected 
physical behavior may be then attributed to the degree 
of localization of the 4f-electrons, i.e. the positions of the 
4f-density of states (4f-D0S) with respect to the Fermi 
level. The position of the 4f-D0S demonstrates the de- 
gree of hybridization between localized 4f and conduction 
bands. The EFG quantity is extremely sensitive to the 
anisotropic charge distributions of the core electrons^, as 
well as to the aspherical electron density distribution of 
valance electrons^, and as a result to the valance elec- 
tronic structure. The EFG, thereby, can serve as a pow- 
erful gauge for measuring such a degree of localization. 

J. Rusz et alj^ very recently calculated Fermi sur- 
faces of the Celns regardless its antiferromagnetic or- 
dering. Their calculations were performed in the local- 
ized extreme limit within the open core trcatmcn t^^'^^ . 
On the other delocalized extreme limit two individual 
groups^ calculated the EFG at the ^^^In site in this 
compound. The calculations of the former— and later- 
groups were performed, respectively, in the antiferro- 
magnetic and nonmagnetic phases employing a similar 
method of the full-potential linearized augmented plane 



waves (FP-LAPW)i^. 

In this paper, we have examined whether one can, us- 
ing an intermediate way, improve the previous results of 
the localized and delocalized limits. For this purpose, we 
have employed the LDA+U schem o^^d^i^^ and then cal- 
culated the EFG's at both "^/n and "°Ce sites. The 
more advanced method of augmented plane waves plus 
local orbital (APW+lo were used to linearize the en- 
ergies. According to our knowledge, this is a first-report 
of the EFG calculations within the LDA+U scheme 
employing APW+lo method for this compound. Our 
spin-polarized calculations demonstrate non-zero EFG at 
^^"^Ce site in the presence of spin-orbit coupling. The 
EFG's, according to our knowledge, had not been previ- 
ously calculated at the cubic ^^°Ce site in this compound. 
The cubic symmetry of the Ce site explains why the 
EFG's were less significant to be previously reported. We 
have emerged, nevertheless, a new result that the EFG's 
are dominated at the ^'^^Ce site by 4f-states and not as 
usual by p-states. The goal of this work is to illustrate 
an approximately linear relationship between the values of 
EFG and density of states (DOS) at Fermi level (Ep), 
viz. EFG (X DOS{Ef). According to our knowledge, 
such a linear relationship had not been previously ob- 
served. We also aim to justify about the tendency of the 
4f-electrons in the ground state of the antiferromagnetic 
Celn^ compound to show their degree of localization. 
We have also found that correlations among 4f-electrons 
infiuence semicore states of 5p-Ce. 



II. DETAILS OF THE CALCULATIONS 

All the calculations in this work were performed in the 
frame work of the density functional theor y ^^d^ (DFT). 
We have taken the generalized gradient approximation^^ 
(GGA) into account for the exchange-correlation func- 



tional. We have employed the full-potential augmented 
plane waves plus local orbital (APW+lo) metho d^^i^^ as 
embodied in the WIEN2k code^i. The mufhn-tin radii 
(Rmt) were chosen to be 2.2 A and 2.8 A for the In and 
Ce atoms, respectively. It has been allowed to be the 
4s 4p 4f orbitals of the Ce and 5d 6s orbitals of In in 
the valance states. The expansion of the wave functions 
inside the spheres in lattice harmonics and in the inter- 
stitial region in plane waves were cut off by the maximum 
eigenvalue of Imax = 10 and the RivirKmax = 7, respec- 
tively. The cutoff for the Fourier expansion of the charge 
density and potential was taken to be Gmax = IGy/Ry. 
We used a mixing parameter of 0.001 in the Broyden's 
scheme to reduce the probability of occurrence of the 
spurious ghostbands. A mesh of 165 special k points 
was taken in the irreducible wedge of the first Brillioun 
zone, which corresponds to the grids of 18 x 18 x 18 
in the scheme of Monkhorst-Pack^^. In order to per- 
form the calculations nearly in the same accuracy for 
the case of magnetic super cell compared to the non- 
magnetic unit cell, the mesh of k points was reduced to 
121 corresponding to 9 x 9 x 9 grids. We diagonalized 
spin-orbit coupling (SOC) Hamiltonian in the space of 
scalar relativistio^ eigenstates using a second- variational 
procedure^! imposing (111) direction on the Ce magnetic 
moments within a cut-off energy of 3 Ry. In order to 
take into account strong correlations of 4f Ce states, we 
have used the LDA+U methodi^'i^'i^. In one side, it is 
generally believed that, the on-site Coulomb repulsion U 
integral, as an input parameter for the LDA4-U calcula- 
tions, depends on the system under study. This means 
that the U parameter can vary for an atom, i.e. Ce in 
this work, from one case to the others. Here, we have 
used a value of 6.2 eV for the U parameter of Ce in 
Celns. Fortunately, this value, U = 6.2 eV, had been 
before calculated^^ for our own case of Celn^. Further- 
more, on the other hand, the value obtained for the U 
parameter depends on the method of calculations. Unfor- 
tunately, the calculations of the U value were performed 
within the linear- mufhn-tin-orbital (LMTO) method^^. 
Although the method of our calculations, i.e. APW-l-lo, 
is different from the LMTO, but we have noticed that for 
another case of face-centered-cubic 7 — Ce the value of 
6.1 eV could give satisfactory results in agreement with 
experiments^. The U=6.1 eV value of fcc-7 — Ce is not 
so far from the U=6.2 eV value of the Ceinz compound. 
Therefore, we did neither more elaboration to calculate 
the Coulomb repulsion U within our APW+lo method, 
nor change this value to optimize it. Another input pa- 
rameter for the LDA+U calculation is the exchange inte- 
gral J. Here, we have used 0.7 eV for the J value. This J 
parameter can be derived using the values of 8.34 eV, 5.57 
eV, and 4.12 eV for the F^, F"*, Slater integrals^^, re- 
spectively, within the following expression, which is valid 
for the f electrons: 
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FIG. 1: "(color online)" (a) Chemical unit cell of Celn^ in 
AuCuz prototype having the lattice parameter of a=4.69 A. 
(b) Constructed conventional supercell with the symmetry of 
face-centered cubic {Pm3m space group) from the chemical 
unit cell by a factor of 2 as a number of chemical unit cells 
drawn in three directions of xyz Cartesian coordinates, (c) 
Primitive rhombohedral unit cell of the constructed conven- 
tional supercell. (d) Magnetic supercell imposing spin order- 
ing of (ID to the Ce moments along (111) axis. 



Finally, since the factor of (U-J) appears in the total 
energy of the LDA+U method instead of U and J in- 
dividually, we set J to zero, and let U be equal to the 
effective value^l of Ue/f = U - J = 5.5 eV. 

III. CHEMICAL, MAGNETIC AND 
ELECTRONIC STRUCTURES 

A. Chemical Structure 

As shown in Fig.[T]Ja), Celns crystalhzes in the space 
group of Pm3m with the binary- fee prototype of AuCu^. 
It is a cubic unit cell, where Ce atoms are located on the 
corners, and the In atoms on the middle of the surfaces. 
The point groups of Ce and In atoms are the cubic m-3m, 
and the non-cubic 4/mmm, respectively. The lattice pa- 
rameter of the Celn^ were measured^^ to be 4.69 A. We 
have used this chemical structure to simulate the non- 
magnetic and ferromagnetic phases. 



B. Magnetic Structure 

It has been experimentally observed that the cerium 
moments in this compound are aligned antiferromagnet- 
ically in adjacent (111) magnetic planes^. Therefore, in 
order to simulate the antiferromagnetic situation, first 
the sides of the non-magnetic unit cell were doubled as 
depicted in Fig. [Ijb) in all the 3 Cartesian xyz direc- 
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tions. We have preserved the space group of PmSm for 
the new magnetic supercell. The preserved fee symmetry 
of the magnetic supercell causes to be reduced the num- 
ber of atoms in the primitive rhombohcdral magnetic unit 
cell, see Fig.[TJc), compared to the conventional unit cell 
shown in Fig.[ljb). Second, we have imposed the antifer- 
romagnetic ordering on the magnetic moments of ^^'-'Ce 
sites. The direction of the 4f spin in the ground state 
electron configuration of Ce atom, i.e. [Xe].4:fh5d^.6s^^, 
were exchanged, i.e. [Xe]AfK5d^.6s^^, alternatively 
with the ordering of (ti) in the (111) direction as shown 
in Fig.md). 



C. Electronic Structure 

We have calculated the density of states (DOS) in the 
lack of both spin-polarization (SP) and spin orbit cou- 
pling (SOC). This calculation has been performed using 
the chemical structure; Fig.[lja). This is what we call it 
nonmagnetic (NM) phase from later on. The SOC was 
then included in two individual steps. The calculations in 
the first step were performed including spin-orbit inter- 
actions among only Ce-clcctrons. Wc refer to the results 
of this calculation by the name of "NM-HSOC(only Ce)" . 
The SOC was then, as the second step, included among 
the In-electrons as well. We call it "NM-t-SOC" phase 
form now on. For sure, the SOC is included in both Ce 
and In atoms in the "NM-I-SOC" phase. Our calculated 
DOS's are in agreement with the previous calculations^, 
and we avoid repeating them here. The SOC not only 
influences the semicore states of 5P Ce and 4d In, but 
also changes the density of states at the Fermi level 
{DOS{Ef))- The SOC splitting in the valance states 
is much smaller than the SOC splitting in the semi core 
states. However, the effects of the SOC would not be ne- 
glected on the results due to the fact that many physical 
quantities, e.g. EFG, are very sensitive to the value of 
the DOS{Ef)- The later point seems to be more signifi- 
cant for the Ce based compounds. A typical LDA/GGA 
calculation produces a sharp density of states for the 4f 
Ce, and situates it at the Fermi level. Thereby, any small 
changes in the sharply located 4f Ce DOS at the Fermi 
level can change more significantly the DOS{Ep), and 
consequently the physical results. 

The DOS's were also calculated for the ferromag- 
netic phase in the lack of SOC, which we call it "FM" 
phase. The result shows that imposing spin polarization 
(SP) causes to be reduced the DOS{Ef) from 141.82 
states/ (spin. Ry) in NM phase to 94.46 states/ (spin. Ry) 
for spin up and to 14.90 states/ (spin. Ry) for spin down 
in the FM phase. Wc have added up and down DOS's 
of the FM phase. Therefore, a total reduction of 32.46 
states/Ry in DOS{Ep) is occurred in going from NM to 
FM phase. Such a reduction affects the physical quanti- 
ties that we shall discuss them in subsequent sections. Wc 
have then included the SOC in the FM phase, which we 
call it "FM-fSOC" phase. The rcsuh shows that includ- 



ing SOC gives rise to a reduction of 14.60 states/Ry from 
FM to FM-HSOC. The DOS[Ef) shows a total reduc- 
tion of 47.06 states/Ry in going from NM to FM+SOC. 
Hence it seems that the SP and SOC can influence phys- 
ical quantities in a similar direction in this compoimd. 

The DOS's were calculated for the AFM state in the 
lack of SOC. The later state is called " AFM(OOl)" phase 
in this paper. In the absence of the SOC there is no pre- 
ferred spatial direction at the cubic ^^°Ce site. In this 
case all the planes, for example (111), (001) and so on, 
are identical. Thus, for simplicity, the antifcrromagnctic 
ordering, (ti), were aligned along (001) axis. The SOC 
interactions were then included in the AFM(OOl) phase. 
One could introduce, at this step, a preferable direction 
in the presence of SOC. However, we would postpone 
setting up the correct direction of the cerium moments. 
The cerium moments were then kept still along (001) 
axis, which we call it " AFM(001)+SOC" phase. The 
cerium moments in the AFM-l-SOC(OOl) phase were not 
changed to the more natural (111) direction, in order to 
avoid mixing the effects of SOC with the effects of spin 
directions in the AFM phase. The rcsuh of AFM(OOl) 
or AFM(001)-hSOC shows that the picks of the semi- 
core Ce 5p and 4d In DOS's were nearly doubled in the 
AFM compared to the NM and FM phases. This is in 
consistent with the fact that the number of electrons are 
doubled in going from the chemical cell to the magnetic 
supercell. Wc have then directed the 4f Ce spins from 
(001) to (111) direction, and recalculated the density 
of states self-consistently including spin-orbit coupling. 
The DOS's qualitatively were similar to those for the 
last case of (001) direction. However, quantitatively only 
changing the direction changes the values of DOS{Ep) 
from 167.47 sates/ (spin.Ry) to 127.98 sates/ (spin.Ry) for 
the spin up, and from 167.24 sates/ (Ry.spin) to 128.24 
states/(spin.Ry) for spin down. Such changes can affect 
the sensitive physical quantities. 

We have included an appropriate correlation among 4f 
Ce electrons with employing LDA+U calculations. The 
4f spin orientations were preserved antiferromagnetically 
along (111) direction, and spin-orbit interactions were 
included in the LDA+U calculations. This constitutes 
our "AFM(lll)-hSOC+LDA+U" phase. Total DOS's 
of the AFM(lll)-hS0C-hLDA4-U phase are shown in 
Figs. [2] The result shows a significant reduction of the 
DOS{Ep) to 28.11 states/(spin.Ry) for both up and 
down spins. Total DOS's are separately illustrated in 
semicore. Fig. [2Jb), and valance. Fig. [2I^c), regions. The 
DOS of the semicore region shows further splitting in one 
of the branches of 5p Cc DOS. This occurs in compari- 
son with the phase of AFM(lll)-hSOC. Our calculated 
total DOS of the later phase, which is not shown here, 
is in complete accord with previous calculations^. One 
may confirm the above mentioned further splitting by 
comparing Fig. 1 given in Ref. d with our Fig. [SJa) or 
(b). One would more clearly observe the further split- 
ting in Fig. m^d). The semicore region in Fig. [^d) is 
restricted to an appropriate interval, where the split- 
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FIG. 2: "(color online)" Total density of states of the AFM(lll)+SOC+LDA+U phase presenting (a) semicore and valance 
states, b) semicore states, (c) valance states, (d) selected interval of semicore states exhibiting further splitting in 5p Ce DOS 
from LDA+U interaction. Fermi level is shown using dashed lines 



ting is occurred. This shows that correlations among 
4f electrons influence not only 4f Ce states directly, but 
also 5p Ce semicore states indirectly in this compound. 
There are two nonequivalent Ce atoms in the magnetic 
cell with opposite 4f spin directions, which are indexed 
as Cel and Ce2. The 4f DOS's of Cel and Ce2 are 
shown in Figs.Ofor both phases of AFM(111)+S0C and 
AFM(111)+S0C+LDA+U. The up and down 4f-D0S's 
of Cel are asymmetric with respect to each other, which 
is the case for Ce2 as well. This is in the case that in 
overall Ce 4f DOS is entirely symmetric regardless in- 
dexes 1 and 2. Therefore, in spite of the fact that Cel 
and Ce2 can individually impose magnetic moment, no 
net magnetic moment is imposed in the whole of the su- 
percell. The 4f DOS of the AFM(111)+S0C phase is also 
shown in Fig. [3] (a). In this case one can more easily com- 
pare it with the later phase of AFM(111)-HS0C-HLDA+U 
shown in Fig. [3] (b). The result shows that the signifi- 



cant reduction in the DOS{Ep) originates mainly from 
the splitting between occupied and unoccupied 4f bands; 
see Fig. [3] (b) . Such a splitting within LDA-I-U treatment 
is due to the effective Coulomb repulsion Hubbard U pa- 
rameter. The 4f Ce DOS piles up in the vicinity of Fermi 
level within our GGA calculations. The later point is 
what one can clearly see for the AFM(111)-|-S0C phase 
in Fig.[3ja). This is in the case that the LDA-I-U calcula- 
tions give rise to be constituted the 4f Ce DOS's far from 
Fermi level. As shown in Fig.[2^b), including repulsion U 
potential energy causes to be shifted down occupied and 
shifted up unoccupied 4f DOS's with respect to Fermi 
level. The separation energy between upper and lower 
4f bands is calculated to be 0.43 Ry. Thus contributions 
of 4f Ce to the value of DOS{Ep) arc highly decreased. 
The receding of the 4f Cc DOS, from the vicinity of Fermi 
level, indicates that the 4f electrons are going to lose their 
itinerant character, and as a result to gain their localized 
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FIG. 3: "(color online)" 4f-DOS of the (a) AFM(111)+S0C 
phase and (b) AFM(111)+S0C+LDA+U phase. Solid (dot- 
ted) curves show up and down DOS's of the Gel (Ge2), where 
label 1 (2) refers to the attributed t (J.) ordering. Fermi levels 
are shown using dashed lines for each of cases. 



character. In order to justify about the tendency of the 
4f-electrons to clarify the degree of their locahzation, one 
can compare Fig. ^c) with Fig. ^h). The comparison 
shows that the 4f Ce states play an important role and 
must not be ignored. The later point may be deduced 
from two subsequent facts. First in energy space, the 
4f states as shown in Fig. ^h) are distributed over an 
energy interval for which other conduction bands, e.g. 
In-states, as shown in Fig. [2ljc) are distributed as well. 
Second in real space, each In atom is surrounded by 4 
Ce atoms as one can see in Fig. [Ija). Consequently, the 
4f states are well hybridized with the other conduction 
bands even after including LDA+U. Thus, the localiza- 
tion is reduced but not vanished by LDA-I-U. Therefore, 
we conclude that not only band-like treatment, but also 
open- core treatment cannot provide satisfactory results 
for the case of localized 4f states in this CelnS compound. 
The former treatment would not be well trusted because 



of putting 4f states right at the Fermi level. Thus the 
band-like treatment overestimates hybridization of the 
4f Ce states with the other valance states. The later 
treatment would not be also well trusted because of con- 
fining 4f states in the core region. Thus the open-core 
treatment usually underestimates hybridization of 4f Ce 
states with the other valance states. 



IV. ELECTRONIC SPECIFIC HEAT 

In this section, to fix ideas it seems advisable to com- 
pare the behavior of the electronic specific heats with the 
behavior of the DOS{Ep) through introduced phases in 
the preceding sec. IIII CI The findings of this comparison 
in the next sec. [B will be of relevance to the goal of this 
paper to realize whether or not such a comparison can be 
generalized to the EFG quantity as well. Therefore, we 
have plotted in Fig. ID^a) the calculated total and 4f Ce 
DOS' s{Ep) versus the discussed phases. The number of 
atoms in all the AFM phases is two times the number 
of other phases. The later point is also the case for the 
DOS's{Ef). We have then divided the DOS's{Ef) by 
2 for all the AFM phases. In this case, one can compare 
the AFM-DOS's{Ef) with the DOS's{Ef) of the other 
phases regardless the number of atoms in the magnetic 
and chemical cells. One performing such a comparison 
can focus only on the magnetic ordering effects through 
all the defined various phases. The behavior of this curve 
constitutes the backbone of this paper. 

We have calculated the electronic specific heats, Cv, in 
the lake of both phonon-phonon and electron-phonon in- 
teractions, and plotted the Sommerfeld linear coefficient, 
7 = Cv/T, in Fig. [Ifb) versus all the phases. The re- 
sult nicely represent the behavior of the total DOS{Ef) 
shown in Fig. [4ja), provided that the calculated specific 
heats per cell were also divided by 2 for the AFM phases. 
One would not be surprised observing such a perfect con- 
sistency between in one side the behavior of the specific 
heats and on the other hand the behavior of DOS' s{Ef) 
going through our defined phases. One would not be so, 
because the specific heats calculations were performed 
taking only electron-electron interactions into account. 
One can analytically omitting phonon interactions eas- 
ily prove^a the formula of Cy /T = l/Sn^klDOSiEF). 
This formula demonstrates a linear relation between 7 
and DOS{Ef). According to our knowledge, it had 
not been analytically established such a relation between 
electric filed gradients (EFG's) and DOS{Ef). There- 
fore, the above sketched strategy will numerically make 
an opportunity in the next section trying to demon- 
strate an approximately linear relation between EFG and 
DOS{Ef). We close this section by reporting the value of 
9.74 mJ / {mol.cell.K^) within our LDA+U calculations. 
This calculated 7 value is almost one order of magni- 
tude less than the experimentally measured^ value of 
130mJoule/(mol.cen.K^). The discrepancy is m agree- 
ment with all the other ab initio calculations for other 
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tensor has been calculated using the following formula—: 
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FIG. 4: "(color online)" (a) Total and 4f-Ce DOS's versus 
a variety of phases discussed in the text at their respective 
Fermi levels. The values of DOS(Ef) are divided by 2 for all 
the AFM phases, (b) Sommerfeld linear coefficient 7 — Cv /T 
in mJoule.mol~^ .cell~^ .K~^ of the electronic specific heat 
shown by crossed-squares versus all the phases. For compar- 
ison the 7 coefficients, as shown by filled circles, are divided 
by 2 for the AFM phases. 



cases in the lack of phonon interactions^ii^i^^. 



V. ELECTRIC FIELD GRADIENT 

The electric field gradient (EFG) is a tensor of rank 2. 
The EFG tensor has only 2 independent components in 
the principle axes system (PAS). The axes of the system 
were chosen such that \Vzz\ > \Vyy\ > \Vxx\- One can 
then only evaluate the main Vzz component of the EFG 

"!7^'"' to determine 



and the asymmetry parameter 77 = 
the two independent components of the EFG in the PAS. 
In this paper we only focus on the Vzz as our calculated 
electric field gradients, since the asymmetry parameters 
are zero for our case. The main component of the EFG 



where radial potential coefficient. V20, within LAPW 
method has been calculated as foUow o^i^'^ : 



V2o{r = 0) 




, Rmt , 

+ 4^^1/(A-)j2(Xi?MT)i1io(i^)- (3) 

K 

The integral yields the EFG contribution of the elec- 
trons inside and over the surface of the muffin-thin sphere 
with radius of Rmt- The summation yields the EFG con- 
tribution of the electrons entirely outside of the spheres. 
The contribution of the electrons inside the sphere is 
called the valance EFG, which we here denote it by V^™'. 
The contribution of the electrons over the surface and 
outside of the spheres is called the lattice EFG, which we 
here denote it by V^?"*. 

Electric field gradients, Vzz's, and their respective 
valance, Vj,™''s, and lattice, Vzz^^s, components were cal- 
culated for all the introduced phases in sec. IIII Gl Our 
results are compared with experiment and other the- 
oretical calculations in Tab. |T1 The anisotropy func- 
tions of Ap{Ep) and Ad{EF) were also calculated and 
listed in Tab. |TT] for all the phases. The result, in 
Tab. shows that our calculated EFG within the phase of 
AFM(111)+LDA + U is in better agreement with experi- 
ment than previous calculations^^. The better agree- 
ment confirms that we have taken more properly correla- 
tions among 4f electrons into account within our LDA-I-U 
calculations. One observes, from Tab. |TT1 that at In site 
the absolute value of Ap{Ef) is one order of magnitude 
greater than the absolute value of Ad{Ep), which is not 
the case at Ce site. For the later case, the absolute value 
of Ap{Ef) is one order of magnitude smaller than the 
absolute value of AdiEp), i.e. \Ap{Ef)\ < \Ad{EF)\. 
We will back to this point soon. What important for 
us here is to study the variation of the EFG and ap- 
propriate anisotropy functions versus discussed phases to 
obtain a relation between EFG and DOS{Ef). There- 
fore, we perform more comparisons throughout Tab. U 
and Tab. |TT] illustratively in Figs. [S] We have plotted 
Vzz in Fig. [5fa) and Ap in Fig. [^b), both evaluated at 
^^^In site, versus all the phases. It is generally believed 
that, contributions to the EFG originating from p states 
dominate^. Consequently, one expects to find a linear 
relation between EFG and Ap, i.e. EFG cx Ap. This 
is what one expects to observe looking at Fig. [3a) and 
Fig. [St^b). It is hard to realize, however, the linear rela- 
tion of EFG (X Ap by comparing these two figures with 
each other. Therefore, in order to exhibit such a linear 
relationship, we have separated 3 nonmagnetic phases 
from the other 6 magnetic phases. The Vzz and Ap were 
respectively shown in Fig. O^c) and Fig. E^d) versus 3 
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TABLE I: The main component of the EFG, Vzz, and its decomposition to valance, 14T , and lattice, components given 

in the units of lO^^V/m^ at both In and ^*°Ce sites for aU the phases discussed in the text together with the calculated results 
within LAPW method by the others as well as experimental EFG. 
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LDA+U 
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In 


Ce 


In 


Cc 


In 


Ce 


NM 


No 


No 




13.209 





13.263 





-0.054 





NM 


Only Ce 


No 




12.466 





12.515 





-0.049 





NM 


Yes 


No 




12.847 





12.538 





-0.051 





FM 


No 


No 




13.018 





13.024 





-0.006 





FM 


Yes 


No 


(001) 


12.843 


0.082 


12.867 


0.105 


-0.024 


-0.023 


AFM 


No 


No 


(001) 


13.027 


0.089 


13.032 


0.086 


-0.005 


-0.003 


AFM 


Yes 


No 


(001) 


12.966 


0.086 


12.968 


0.066 


-0.002 


-0.002 


AFM 


Yes 


No 


(111) 


12.857 


0.029 


12.961 


0.040 


-0.004 


-0.011 


AFM 


Yes 


Yes 


(111) 


12.431 


-2.863 


12.442 


-2.889 


-0.011 


0.026 


®AFM 


Yes 


No 


(111) 


12.490 




12.540 




-0.050 










unknown 


11.6 













®Ref.3_ 
^Ref.X 



TABLE II: Valance p and d anisotropy functions, Ap{Ef) and Ad{EF), evaluated at In and ^''"Ce sites for the discussed phases 
in the text. 



Phase 




Inte. 








Ap 






Ad 




SOC 




LDA+U 




In 




Ce 


In 




Cc 


NM 


No 




No 




0.1066 







-0.0060 







NM 


Only Ce 




No 




0.1040 







-0.0061 







NM 


Yes 




No 




0.1043 







-0.0061 







FM 


No 




No 




0.0373 







-0.0050 







FM 


Yes 




No 


(001) 


0.0369 







-0.0052 







AFM 


No 




No 


(001) 


0.0373 




0.00015 


0.0025 




0.00130 


AFM 


Yes 




No 


(001) 


0.0372 




0.00005 


0.0024 




0.00130 


AFM 


Yes 




No 


(111) 


0.0370 




-0.00010 


0.0024 




-0.00085 


AFM 


Yes 




Yes 


(111) 


0.0361 




0.00715 


0.0026 




0.01880 



nonmagnetic phases. Similar results, Vzz and Ap, for 
the magnetic phases were separately shown in Fig. [S{e) 
and Fig.[5l^f). Now one can among non magnetic phases 
clearly observe the linear relation between EFG and Ap 
by comparing the behavior of Vzz, Fig. HJc), with the 
behavior of Ap shown in Fig. [Sl^d). They behave simi- 
lar to each other trough non magnetic phases. This is 
also the case for the magnetic phases as well, since the 
Vzz shown in Fig. [S{e) behaves similar to Ap shown in 
Fig. m^f) through the magnetic phases. Now time seems 
apt to intuit that there is an approximately linear relation 
between EFG and DOS{Ef), i.e. EFG oc DOS{Ef). 
The relation can be realized, if the behavior of the to- 
tal DOS{Ep) shown in Fig. [HJa) is compared with the 
behavior of the Vzz shown in Fig. \5lc) and Fig. [Sje). 



The comparison yields the result that the EFG is ap- 
proximately proportional to the DOS{Ef)- The result is 
obtained, because the Vzz and the DOS{Ep) vary simi- 
larly, going through all the phases. This is analogous to 
the used strategy in sec. IIVI to realize the linear relation 
of Gv/T = l/3Tr^klD0S{EF). Our resuh shows that 
the proportional constant is negative for the non mag- 
netic phases, while it is positive for the magnetic phases. 
The former constant is negative, because the Vzz and 
DOS{Ef) inversely vary through non magnetic phases. 
They, however, vary directly versus magnetic phases re- 
sulting to the positive constant. 

Furthermore, one can also conclude that the value 
of 4f-D0S(EF) can significantly influence the value of 
EFG. One may confirm this conclusion due to a simi- 
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(C) 



(d) 



13.0- 
12.8- 
12.6- 

12.4- 

0.0372 
0.0368 
0.0364 
0.0360 



'■■»■■' 



(f) 



c9 

^ ^ <r 'T 
/ / / 



FIG. 5: "(color online)" (a) The calculated EFG at the ^^^In site and (b) the In 5p anisotropy function, Ap{Ef), for all the 
nonmagnetic and magnetic phases discussed in the text, which they are also represented in (c) and (d) for only nonmagnetic 
phases as well as in (e) and (f) for only magnetic phases. 



larity between the behavior of 4f DOS{Ep) and total 
DOS{Ep) shown in Fig. 2] (a). The foregoing result of 
EFG cx DOS{Ep) ensures that the EFG can be also pro- 
portional to the 4fD0S{EF); i.e. EFG cx DOS'^f{EF). 

The later proportionality between EFG and 4f 
DOS{Ef) makes more crucial the method of treatment 
with 4f-Ce electrons. The later point describes why the 
value of EFG at ^^^In site within the LDA-f U calculation 
is smaller than those obtained within all the other cal- 
culations performed in the lack of LDA+U interactions. 
The description can be provided taking the splitting of 
0.43 Ry shown in Fig. [3[b) into account between occu- 
pied and unoccupied bands due to the LDA-I-U treat- 
ment. The splitting gives rise to be shifted downwards 
the 4f DOS from the vicinity of Fermi level, and as a 
result to be reduced the value of 4f DOS{Ef)- The 
later reduction together with the discussed relation of 
EFG oc DOS^^ {Ef) provides the satisfactory descrip- 
tion why the EFG is reduced within LDA-I-U calcula- 
tions. All these support our previously concluded result 
in sec. IIII C1 concerning the important role of 4f electrons 
in this compound. 

At Ce site due to its cubic point group the EFG's 
are zero, as listed in Tab. [H for all non magnetic phases 
and FM phase. For FM-f-SOC phase and all the AFM 
phases, however, our result shows non zero EFG values 



at this site. The non zero values for the EFG's origi- 
nate from the fact that SOC or magnetic ordering can 
give rise to a little bit deviation from cubic symmetry. 
One expects that the small deviation gives rise to small 
EFG's. Our result, in Tab. HI confirms the later point 
apart from the last phase of AFM(lll)-t-SOC-f LDA+U. 
For the later AFM(lll)-hSOC-hLDA-hU phase the EFG 
at Ce site is also practically small, but it is 2 order of 
magnitude larger than the other phases. To find the 
source of such a discrepancy, we follow our last strategy 
looking at the behavior of EFG and anisotropy functions 
through corresponded phases. The calculated EFG at 
Ce site can be compared with the Ap versus magnetic 
phases (apart from FM phase for which EFG is exactly 
zero) using Fig. [HI (a) and (b). The comparison does not 
show similar behavior for the EFG and Ap along the 
magnetic phases. To find the reason why they do not 
show similar behavior, now we come back to the men- 
tioned point that \/S.d{EF)\ > \Ap{Ef)\ at the Ce site; 
see Tab. [ITl One first suspects, due to the larger value 
of \Ad{EF)\ than \Ap{Ef)\, that the behavior of EFG 
might be similar to the behavior of Ad{EF)- The behav- 
iors of EFG and Ad{EF) are compared in Fig. [5] (a) and 
(c). However, they also do not show similar behavior. 
Therefore, we could not reproduce the behavior of the 
EFG using either Ap{Ef) or Ac?(£'i?). For more rcaliza- 
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FIG. 6; "(color online)" (a) The calculated EFG at the 
^''"Ce site, (b) the Ce p anisotropy function, Ap{Ef), and 
(c) the Ce d anisotropy function, Ad{EF), for the phases dis- 
cussed in the text for which nonzero EFG's have been cal- 
culated at the ^*°Ce site, viz. 
AFM(lll)-hSOC-hLDA-^U. 



ranging from FM-fSOC to 



tion, wc have decomposed the EFG at Ce site into its 
valance contributions. The results are shown in Figs. [7| 
for spin up and spin down. In the Figs. [71 we have inset 
a figure to show sum of up and down spins for each of 
the valance contributions individually. The result, com- 
paring Fig. [7] (a) and (b), shows that contributions of 
d-states to the EFG is one order of magnitude less than 
contributions of p-states. This is in the case that, as 
expressed before, \Ad{Ep)\ is one order of magnitude 
larger than |A|3(£'j?)|. Therefore, even for the case of 
|Ad(£'i?)| > \Ap{Ep)\, contributions to the EFG origi- 
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FIG. 7: "(color online)" (a) Decomposed up, t, down, 
I, and up plus down, t + |(see insets), valance contri- 
butions of EFG evaluated at ^'^°Ce site for the magnetic 
phases discussed in the text ranging from FM-I-SOC to 
AFM(lll)-(-SOC-hLDA+U. (NB: There are no any EFG con- 
tributions at Ce site for FM phase, in the lack of SOC, due 
its cubic symmetry.) 



10 



nating from p-states compared to the d-states dominate. 
This result is in c omp lete accord with the calculations 
performed in Ref. |3J for transition metals. However, a 
different and important result can be emerged from com- 
paring Fig. [7] (c) with Fig. [6] (a). The result is this that 
the EFG mainly originates from f-states at Ce site. We 
have added all the valance contributions to the EFG in 
Fig. [3 (d). The result shows that the behavior of EFG 
shown in Fig. [6] (a) is similar to total valance EFG con- 
tributions shown in Fig. [7| (c). Lattice contributions to 
the EFG also cannot change the similarity, because they 
are negligible as listed in Tab. [H Therefore contributions 
to the EFG originating from f states dominate at Ce site. 
For sure in Fig. [7] (e) we have subtracted contributions 
to the EFG of f-states from total valance EFG contribu- 
tions. The result, Fig. [7] (e), represents the behavior of 
p-anisotropy function shown in Fig. [6] (b) . Such a repre- 
sentation reconfirms our last conclusion. 



VI. CONCLUSION 

We have investigated the variations of electric field gra- 
dient (EFG's) and their valance contributions as well as 
their anisotropy functions at both In and Ce sites in the 
Celn^ for a variety of circumstances. For each of the cir- 
cumstances, density of states at Fermi level {DOS{Ep)) 
is calculated. We have found that comparing the behav- 
ior of the EFG with the behavior of DOS{Ep) versus the 



applied circumstances may be a useful strategy to emerge 
physical properties. We have found within such a strat- 
egy that the main component of the electric field gradi- 
ent, Vzz^ is approximately proportional to the value of to- 
tal density of states at Fermi level {DOS{Ep)) as well as 
4/ — DOS{Ef). Despite anisotropy function of d-states 
is larger than the one of p-states at Ce site, contributions 
to the EFG originating from p-states dominate compared 
to d-states. This is in the case that, however, the EFG's 
are dominated at the ^^'^Ce site by the 4f-states, while 
they are as usual dominated by p-states at the In site. 
The result shows that 4f Ce states are hybridized with 
conduction bands and play an important role which must 
not be ignored. Thus we have within our electronic struc- 
ture calculations predicted that neither band-like treat- 
ment, nor open-core treatment can provide satisfactory 
results for this case. Our LDA+U calculations arc in bet- 
ter agreement with experiment indicating the fact that 
correlations among 4f electrons are taken more properly 
into account. The correlations among 4f electrons in- 
fluence not only 4f Ce states directly, but also 5p Ce 
scmicore states indirectly in this compound. 
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